

pacman::p_load(here, fst, data.table, ggplot2)

raw_cems = read.fst(
  here("Data/electricity/raw_cems_2017.fst"),
  as.data.table = TRUE
)

raw_cems[,.(tot_gload = sum(gload, na.rm =T)),by = .(op_hour)] |> 
  ggplot(aes(x = op_hour, y = tot_gload)) + 
  geom_line()
         
